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We show that a self-interacting neutrino gas can spontaneously acquire a non-stationary pulsating 
component in its flavor content, with a frequency that can exactly cancel the “multi-angle” refractive 
effects of dense matter. This can then enable homogeneous and inhomogeneous flavor conversion 
instabilities to exist even at large neutrino and matter densities, where the system would have 
been stable if the evolution were strictly stationary. Large flavor conversions, especially close to a 
supernova core, are possible via this novel mechanism. This may have important consequences for 
the explosion dynamics, nucleosynthesis, as well as for neutrino observations of supernovae. 
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I. INTRODUCTION 

Inside a supernova (SN), neutrino densities are so high 
that neutrino flavor oscillations are affected not only by 
ordinary matter, but also by the neutrinos themselves 
[4]. Neutrino-neutrino interactions lead to highly corre¬ 
lated collective flavor conversions and unexpected effects, 
which completely change the physics of supernova neu¬ 
trinos IM23]. See refs. mils] for recent reviews. In this 
work, we show that time-dependent fluctuations lead to 
a novel effect that can enable flavor conversion deeper in 
the SN than previously realized. 

Above the neutrinosphere, in the absence of collisions, 
the dynamics of a dense neutrino gas is characterized in 
terms of “density matrices” in flavor space, p(t, x, E, v) 
for neutrinos with energy E and velocity v at position x 
and time t. These obey the kinetic equations [21] |26l |27] , 

i{dt+-v ■'V)q= [H,p] . (1) 

On the left-hand-side (l.h.s.), the first term accounts 
for explicit time dependence, while the second, propor¬ 
tional to the neutrino velocity v, is the drift-term due to 
neutrino free-streaming. On the right-hand-side (r.h.s.), 
H(t,x, T^,v) is the Hamiltonian matrix in flavor space 
containing the neutrino mass-square matrix and poten¬ 
tials due to matter and neutrinos. 

Flavor evolution of the dense neutrino gas, as governed 
by eq. has a highly complex structure. It depends on 
the 4 time and space coordinates, the 4 energy and ve¬ 
locity coordinates (with |v| = 1, in our ultra-relativistic 
approximation), as well as the flavor states of all neutri¬ 
nos. In order to reduce this complexity, symmetries in 
the neutrino flavor evolution have often been assumed. 
For neutrinos in a SN environment, all of previous lit¬ 
erature is based on the assumption that the evolution 
is stationary, i.e., there is no explicit time-dependence, 
or only a slow/small time-dependence that does not sig¬ 
nificantly affect the flavor evolution. Additionally, under 
the assumption of a spherically symmetric neutrino emis¬ 
sion, the dynamics reduces to a one-dimensional evolu¬ 


tion along the radial coordinate. This is the rationale 
behind the often-used “bulb model” mm- 

These symmetry assumptions, viz., temporal sta- 
tionarity and spatial homogeneity, have been recently 
criticized because self-interacting neutrinos can spon¬ 
taneously break these space-time symmetries. Indeed, 
studies on simple toy-models show that the translation 
symmetries in time [23 [29] and space [3QH3^ are not 
stable. Even tiny inhomogeneities may lead to new fla¬ 
vor instabilities [3Qll3^ that can develop also at large 
neutrino densities, as above the SN core, where oscil¬ 
lations are otherwise expected to be suppressed due to 
synchronization. However, large neutrino densities in 
a SN are typically accompanied by a large matter den¬ 
sity [34| , which produces “multi-angle matter effects” [35] 
that suppress both homogeneous and inhomogeneous in¬ 
stabilities. The current understanding is then that neu¬ 
trinos cannot change their flavor too close to the SN core. 

Flavor conversions at small distances from the SN core 
would have major consequences for SN explosions, nu¬ 
cleosynthesis, as well as neutrino observations of nearby 
SNe. If conversions are possible below the shock radius, 
neutrinos can provide a net positive energy to the shock 
and assist SN explosions [361139] . Similarly, the neutron- 
to-proton ratio can be changed deeper inside a star, af¬ 
fecting the yield of heavier nuclei created through the r- 
process gnun]. Also, in order to interpret any potential 
observation of neutrinos from SNe, current and proposed 
neutrino experiments depend crucially on understanding 
where and how the flavor-dependent neutrino fluxes have 
converted to each other [25l[42H^ . Now that Gd-doping 
in Super-Kamiokande [45] is approved [46] , the imminent 
observation of the diffuse background of SN neutrinos 
may raise this issue m, even without a Galactic SN. 

In the following, using linear stability analysis we show 
the presence of an unstable pulsating mode that leads to 
flavor conversion at high neutrino and matter density. 
The key insight is that the frequency of pulsation can 
undo the phase dispersion due to a large matter den¬ 
sity. As a result, flavor instabilities, which would have 
grown only if matter effects were small, can now develop 
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at large neutrino and matter densities. Then, to demon¬ 
strate that this linear instability survives in the non¬ 
linear regime, we numerically calculate the flavor evo¬ 
lution in a simplified model and show that flavor conver¬ 
sions indeed occur at large neutrino and matter densities 
when there are space and time-dependent fluctuations. 
Finally, we discuss the implications for SN neutrinos and 
conclude. 


II. LINEAR ANALYSIS FOR A GENERAL 
SCENARIO 


Assuming that the neutrinos are initially in flavor 
eigenstates, their density matrices p(t, x, F’, v) can be 
written in a 2 -flavor framework as 


Q = 


Tr(g) , ^ 

2 2 ^ 




( 2 ) 


to linear order in S'(t,x, F,v) [48]. The quantity 
^(t,x, F,v) is the energy and angular distribution of 
neutrinos from the source and is an arbitrary nor¬ 
malization constant, with dimensions of number density, 
for making S dimensionless. A non-zero off-diagonal ele¬ 
ment S represents flavor conversions. For antineutrinos, 
q{E) = —g{—E), extending the physical range of E from 
—oc to -hoc. The Hamiltonian for the flavor evolution is 

H = + ^/2GpMi + y/2Gp J dT (1 — v • v , (3) 

where sans-serif (quantities are 2 x 2 matrices in flavor 
space. Namely, is the neutrino mass-squared matrix, 
while \/2GfNi and the last term on the r.h.s. appear due 
to forward scattering on matter m and neutrinos mia, 
respectively. The integral is over all neutrino energies 
and velocities, i.e., f dV' = dE' E'^ J dw' 

In a non-isotropic neutrino gas, as in the case of neutri¬ 
nos streaming off a SN core, there is a net neutrino cur¬ 
rent so that neutrinos moving in different directions, i.e., 
with different v, acquire different phases via the velocity- 
dependent terms, i.e., (1 — v • v^) in H and v • V from the 
drift-term in eq. 0. These are multi-angle effects, that 
arise due to the current-current nature of the low-energy 
weak interactions and the source geometry. Typically, 
they inhibit the collective behavior of the flavor evolu¬ 
tion, but can also lead to flavor decoherence [50H52]. 

Using eqs. <0 and 0 in eq. 0 > and taking a vanishing 
mixing angle, we find the equation for flavor evolution, 

i[di + V • V)S = [ — cj + A-f/i/ dr'{l — V • S 

- IjJ dT'{l-v ■v')g'S', 

( 4 ) 


where the relevant energy scales are the neutrino oscil¬ 
lation frequency in vacuum uj = Am^/( 2 F), the mat¬ 
ter potential A = y/^GpU^, and the neutrino potential 


li = y/2Gpnj^. Note that /i always appears in product 
with S', making the precise choice of immaterial. 

Let us consider the evolution of S along the radial dis¬ 
tance r, while Fourier decomposing it in t and the spa¬ 
tial coordinates transverse to f, viz., Yj^. We will take 
the spectrum g x, F, v) to be independent of time and 
space, i.e., g = g'{E^v)^ so that it does not 

get Fourier transformed. Explicitly, 

/ -|-oo 

-CXD 

where is the Fourier coefficient of a flavor 

evolution mode with temporal pulsation p and inhomo¬ 
geneity wavevector k . Inserting this ansatz into eq. Q , 
using V • V = and dividing by the radial 

velocity we find an eigenvalue equation for Qp^\^{E, v). 


—uj X — p — ' k 

Vj. 



Q 


p,k 


ii 



- V • v')g'Qp,\^ , 


(6) 


where A = A + /i / dr'{l — v • v')g' encodes “matter” 
effects from both matter and neutrinos. Note that in the 
linear regime, different Fourier modes are not coupled. 
A growing solution to this equation, with Im(U) > 0, 
signals that there is an instability. 

For a stationary system one finds growing solutions 
even at a large neutrino density, if inhomogeneities are 
present, i.e., k 7 ^ 0 , as long as A ^ /i [30l|3ll[34]. In 
a SN however, A is also large when p is large and these 
instabilities are typically not realized [34] . The reason for 
this is clear: One cannot simultaneously obtain growing 
solutions for all velocities, because the inhomogeneous 
term, the matter term, and the /i-dependent neutrino- 
neutrino interaction term on the r.h.s. are all large, i.e., 
-k. A, /i ^ cj, but have different velocity-dependences 
which cannot completely cancel against each other. 

The non-stationary system has an innocuous-looking 
but important difference with respect to the stationary 
system. Non-stationarity lowers A by p, i.e., A ^ A—p, as 
also seen in ref. [29] . More importantly, a fact not realized 
so far, this p 7 ^ 0 term has the same multi-angle depen¬ 
dence as A. Therefore, if one allows for a non-stationary 
solution, the neutrino system with a pulsation p A can 
undo the phase dispersion due to a large matter term 
for all velocities. Thus, one can find growing solutions, 
with lm( 0 ) > 0 , to the eigenvalue equation [eq. (|^], as 
previously for the X p scenario. The eigenvalues are 
identical to those in Sec. 4.2.4 of ref. [34] with the shift 
A ^ A — p. These solutions are highly oscillatory in 
space (k 7 ^ 0 ) and time (p 7 ^ 0 ), and would lead to flavor 
averaging. This is our main result. 

In Fig.[^ we illustrate this general idea in the context 
of SN neutrinos. The thin black line shows a SN den¬ 
sity profile, while the thick colored lines schematically 
show where flavor instabilities with certain pulsation p 
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Figure 1: Schematic of the SN matter potential (thin line) 
and densities where neutrino flavor composition is unstable 
to small oscillations with spatial inhomogeneity wavenumber 
k and frequency of time-dependence p (thick lines) [53]. Time- 
dependent fluctuations can be unstable at higher matter den¬ 
sity and lead to neutrino flavor conversions deeper in a SN. 

and wavenumber /c = |k| can grow, i.e., have Im(f^) > 0. 
The inhomogeneous instability, i.e., with k ^ 0, is al¬ 
ways below the homogeneous k = 0 instability [34|, and 
never occurs for the physically available matter density. 
However, p X can raise these k = 0 and k 0 instabil¬ 
ities, making them unstable at low radii. They can now 
develop at large /i and A, i.e., at a small radius, with a 
temporal oscillation of frequency p. This can happen 
for both normal and inverted neutrino mass ordering. 

How does this linear instability evolve when linear the¬ 
ory is no longer appropriate? What is its impact? To 
answer these questions, we must numerically solve the 
equations of motion (EoMs), eq. in the fully non¬ 
linear regime. This is what we do next. 

III. NON-LINEAR ANALYSIS FOR THE 
TWO-BEAMS MODEL 

The new effect, relevant to SN neutrinos, is that a 
stationary system, which is stable to both homogeneous 
and inhomogeneous perturbations at large fi and A, be¬ 
comes unstable therein when non-stationarity is allowed. 
This requires simulating a system with temporal non- 
stationarity, spatial inhomogeneity, and multi-angle mat¬ 
ter effects, which is extremely challenging and has not 
been attempted so far. Here, we present the first simu¬ 
lation with these three features. 

The simplest model that can accommodate the re¬ 
quired features is the neutrino “line model” |30l [32] |34] . 
In this model, one considers monochromatic neutrinos 
emitted in two directions, “L” and “i?”, from an in¬ 
finite plane at 2: = 0. Assuming translational invari¬ 
ance along the ^-direction, the flavor evolution along 
2; > 0 can be characterized on the two-dimensional plane 
spanned by the x and ^ coordinates. The neutrino emis¬ 
sion modes L and R are labeled in terms of their veloc- 


ities, i.e., 0, = (cosi^^,0,sini?^), where 

G [0,7r] is the emission angle, and similarly for v^. 
Thus, for the L mode, the differential operator on the 
l.h.s. in eq. 0 takes the form 

+ W • V = 9^ + , (7) 

while the Hamiltonian in eq. 0 becomes 

Hl = - 2 + /^(l “ ' ^r) [(1 + ^)Qr - Qr] : 

(8) 

where is the diagonal Pauli matrix, and e is the 
neutrino-antineutrino asymmetry, i.e., 1 + e = (n^^ — 
~ with being a non-electron flavor. 

The normalization = rip^ — rip^ is used to define 
11 = V2Gf -npj. 

The differential operator in eq. 0 shows that the flavor 
evolution is determined by a partial differential equation 
in one temporal and two spatial dimensions. By Fourier 
transforming the EoMs in t and x, as in eq. ([^, one ob¬ 
tains a tower of ordinary differential equations in the 2:- 
coordinate for the different Eourier modes Qp j. with tem¬ 
poral pulsation p and spatial wavenumber k. In the non¬ 
linear regime, the EoMs for the different Eourier modes 
have a convolution term due to interactions between the 
different modes [32|. 

Vz L — Ri assumed in refs. |30l|32], even a very 
large matter term A can be rotated away from the EoMs 
by studying the flavor evolution in a suitable co-rotating 
frame [6|. Conversely, if Vz^l 7^ f^e matter term 
leads to frequencies Xjvz^p fhe L mode and XjVz^R 
for the R mode. Their difference cannot be removed and 
gives multi-angle matter effects [35]. If A ^ /i, this phase 
difference between L and R modes is so large that it 
suppresses the self-induced flavor conversions from both 
k = 0 and k 0 instabilities [SlISS]. Conversely, if one 
allows for a non-stationary solution, the neutrino system 
selects the pulsations p A, that compensate the phase 
dispersion due to a large matter term, and generate grow¬ 
ing instabilities, in particular at small-scales associated 
with spatial inhomogeneities. 

To quantitatively illustrate this claim, we take the 
source at 2; = 0 to emit only and Pg, with an fac¬ 
tor of two excess of r^ over Pg, i.e., e = 1. We choose 
0 = and a normal mass ordering, i.e., u; > 0, 

but the result would be similar for the inverted order¬ 
ing, i.e., uj < 0. The overall frequency-scale is set by 
uj = 1. A large /i = 40 is chosen, so that oscillations 
are suppressed in the homogeneous case, as in a SN. We 
take the L and R modes to have two different angles 
i^R = 57r/18 and so that a large matter 

potential, A = 4 x 10^, suppresses the inhomogeneous 
modes, mimicking the similar effect in a SN. 

In order to make k and p dimensionless, k is expressed 
in multiples of the vacuum oscillation frequency, i.e., 
k = while p is expressed in multiples of the matter 
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Figure 2: Amplitudes of flavor conversion at distance z in our two-beams model, log^gj^^^ (z)|, for the Fourier modes of 
inhomogeneity along the transverse direction. Time-dependent pulsations of frequency p — \ and large neutrino and matter 
densities are included. Inset: Linear growth rates Im(Q)/ct; for a larger range of (note the log-scale). Modes ~ 10^ grow 
first and quickly excite all other modes, leading to large flavor conversion. 


potential A, i.e., p = UpX. For simplicity, we limit our¬ 
selves to the Up = 1 mode. In this way, in the non-linear 
regime we only have to consider the convolution among 
the different Fourier modes associated with spatial inho- 
mogeneities, ~ EjJ((l + e)QR,n,-j, “ QR,n,-jJ,QL,jJ, 
and analogously for the R mode [32]. We include modes 
up to Uf. = 600 but ensure that nj. > 400 remain empty. 
This trick avoids “spectral blocking” that leads to a spu¬ 
rious rise of the Fourier coefficients at large rij. due to 
truncation of the tower of equations [54]. To seed the 
spatial inhomogeneity, we use numerical noise of O(10~^) 
for all modes. 

Fig-i shows the flavor dynamics of this two-beams 
model. Amplitudes of flavor conversion, log^ol^n^ (^)|5 
are shown at distance ^ for various Fourier modes rij^ of 
inhomogeneity along x, in the presence of time-dependent 
fluctuations with frequency p = X. As linear theory pre¬ 
dicts [Im(f2) vs.n/j. shown in the inset], modes around 
Til. 100 are the most unstable and grow first with 
the predicted rate Im(f2)/ci; 3. However, while only 
the modes O(IO^) are unstable according to lin¬ 

ear analysis, a cascade in Fourier space develops to both 
smaller and larger Uf. due to the convolution enforced by 
the neutrino-neutrino interaction. Modes with rij. < 10^ 
also grow fast when nonlinearity sets in. The flavor com¬ 
position begins to oscillate with many frequencies and 
appears to be “averaged out”. 


Thus, we find that neutrinos can change flavor at large 
A and /i, if non-stationary solutions with frequency p X 
are allowed. We expect that the cascade in Fourier space 
leads to “flavor decoherence”, or approximate equilibra¬ 
tion between all flavors [55]. In a SN, this may have 
important consequences, which we discuss below. 


IV. DISCUSSION AND CONCLUSIONS 

An important question is if this effect can be impor¬ 
tant in a SN. Here, we provide a back-of-the-envelope 
estimate. To aid shock revival, flavor instability has to 
occur below the shock-front at r 200 km and above 
the gain radius at r 100 km (see, e.g.. Fig. 4 in 

ref. [38]). One thus needs pulsations of high frequency 
p ^ A ^ 10^~^km~^ ^ 3 X 10^~^^Hz. Between r 
100 km to 150 km, a typical instability with growth rate 
Im(U)/cc; 3 then grows by ^ 60 e-foldings, i.e., a factor 
of 10 , for 15MeV neutrinos with oj ^ 0.4 km , as¬ 
suming constant growth. Can such high-frequency fluc¬ 
tuations occur in a SN with even tiny amplitudes? In 
this context, we find it intriguing that pair-correlations 
of the neutrino field, which are many-body corrections 
to the single-particle density matrices, lead to relative 
number fluctuations of a size ^ {X^/E)‘^ ^ 10~^^, 
where P 10 c is the typical speed of ordinary matter 
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in SN, which oscillate with a frequency ^ 2E r\j lO^^Hz 
for 15MeV neutrinos, as shown in ref. [23]. However, in 
a realistic SN the density behind the shock, though often 
much flatter than shown in Fig.[^ is not exactly constant 
and the growth rates decrease when p ^ X. On the other 
hand, then nearby p modes are excited and nonlinear 
coupling of modes makes all instabilities grow. More de¬ 
tailed studies are needed to study these effects [56|, the 
impact on nucleosynthesis, and overall flavor conversion. 

Potentially, the consequences of our finding augur an¬ 
other paradigm-shift in the understanding of self-induced 
conversions and on their impact on the SN dynamics. 
The possibility of low-radii conversions behind the stalled 
shock-wave during the accretion phase, suppressed by the 
large matter term in the stationary and homogeneous 
case EMU, implies that the flavor dynamics may need 
to be taken into account in the revitalization of the shock- 
wave |36l [38]. Also, the impact on nucleosynthesis in a 
SN would be important IlDlllI]. With flavor equilibra¬ 
tion, the interpretation of observed SN fluxes may also 
become simpler ESI- 

However, the possibility of flavor conversions close to 
the neutrinosphere in a SN also questions the assumption 
that flavor conversions safely occur outside it. This as¬ 
sumption allowed one to replace the full Boltzmann equa¬ 
tions, containing both oscillations and scatterings, with 
the flavor oscillation equations for free-streaming neutri¬ 
nos. In a non-stationary situation, this assumption may 
no longer be guaranteed and imply the necessity to si¬ 
multaneously perform the neutrino transport and flavor 
evolution m- This is a formidable problem that would 


require new computational techniques. Furthermore, it is 
possible that these instabilities (inhomogeneity and non- 
stationarity) appear in a regime where the coarse-grained 
description adopted using density matrices [2] is insuffi¬ 
cient. Although we have used this standard description 
here, as a first step, this is a more fundamental aspect 
that needs further studies. 

In conclusion, we have presented the first study of non¬ 
linear effects of non-stationarity in a dense neutrino gas. 
We have pointed out novel temporal instabilities that can 
dramatically affect flavor evolution, and raise the possi¬ 
bility of self-induced flavor conversions deep in a SN. The 
discovery of the role of symmetry-breaking in the flavor 
evolution of SN neutrinos is opening completely new di¬ 
rections of investigations. We foresee that many surprises 
are still in store. 
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